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Abstract 

The Rabi Hamiltonian describes a single mode of electromagnetic radiation interacting with 
a two-level atom. Using the coupled cluster method, we investigate the time evolution of this 
system from an initially empty field mode and an unexcited atom. We give results for the atomic 
inversion and field occupation, and find that the virtual processes cause the field to be squeezed. 
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I. INTRODUCTION 



The Rabi Hamiltonian plays an important role in quantum optics. It describes a two- 
level atom interacting with a single mode of quantized electromagnetic radiation via a dipole 
interaction [|I| . It also finds wider application, describing a spin interacting with phonons in 
NMR 10], for example. It is also related to the static Lee model in field theory [^]. 

The field mode is described by bosonic annihilation and creation operators, b and 
respectively, which obey the usual commutation relation, 

[b,b^]=l. (1) 

The two-level atom is described by the pseudo-spin operators 

^^^o-J' ^^^(ooj' ^"^(20)^ 

a'^lia^+a-). a' = '- (a' - a*) . (3) 
With these definitions the Rabi Hamiltonian is given by 

H = -uooo' + Lob'^b + g{a+ + a') {b^ + b) . (4) 

There is a conserved parity 11 associated with the Hamiltonian, 

H = exp (zttA^) , N = b^b+l{a' + 1) . (5) 

There is no proof that this Hamiltonian is integrable, although suggestive evidence does 
exist [^]. Consequently, investigation of the Rabi system requires approximations to be 
made. 

The most widely used approach is to make the rotating wave approximation (RWA). 
This was first applied by Jaynes and Cummings [^], yielding the Jaynes-Cummings Model. 
Under the RWA, one neglects the "virtual terms" a^b'^ and cx^b, also known as the counter- 
rotating terms. This leads to the excitation number as well as parity H being conserved, 
and renders the model soluble by a series of 2 x 2 diagonalisations. 

Despite its frequent use in quantum optics the actual validity of the RWA for specific 
applications is usually highly questionable. For example, we know that the energy spec- 
trum of the Rabi Hamiltonian can be approximated by its JC counterpart in the RWA only 
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for sufficiently small values of the coupling strength and that the width of this range 
decreases as one proceeds higher up the spectrum. There is a further concern regarding 
the RWA. Ford and O'Connell p[ investigated the system of a charged harmonic oscillator 
interacting with a reservoir consisting of an infinite number of oscillators within the RWA. 
They showed that the spectrum has no lower bound for all models of physical interest. This 
constitutes a violation of the second law of thermodynamics, as we can take energy from 
the oscillator bath (modeling the environment) without producing an effect upon it, i.e., we 
can still remove an infinite amount of energy from it. 

The energy level spectrum of the full Hamiltonian has been investigated by several au- 
thors. The simplest approach is to use the configuration-interaction (CI) method, equivalent 
to a large-scale diagonalisation in a suitably defined finite subspace of the full Hilbert space. 
This method has been used both by Graham and Hohnerbach [^] and by Kus in in- 
vestigating possible quantum chaos signatures of the model. Lo et al. have given an 



analysis of the validity of the CI method. Reik and others [|T2[ have adapted Judd's method 
|TT| for the Jahn- Teller system for use with the Rabi Hamiltonian. Here, the Hamiltonian 



is translated into the Bargmann representation [13| and solutions of the resulting differen- 
tial equations are sought. Whereas Judd originally used a power series Ansatz, Reik uses 
a Neumann series. For certain couplings, the Neumann series terminates, giving isolated, 
exact solutions known as Juddian solutions. Elsewhere the series gives a useful, convergent 
approximation. The Juddian points are valuable for comparison of approximate techniques. 



Variational results have also been provided by Bishop et al. and by Benivegna and 
Messina |T^. The latter method also permits perturbative corrections, allowing the exact 
results to be approached. 

The time evolution of the Rabi Hamiltonian has been of considerable interest for a long 
time. We mention in particular the pioneering studies of Shirley [|1^. Somewhat later, 
Eberly and co-workers [jl^ were the first to fully demonstrate the rich time evolution of 
this model. They worked within the RWA and discovered the collapse-revival nature of the 



subsequent evolution when the system is started from a coherent field [0. The impact of the 
non-RWA terms on the time evolution has been investigated using various techniques. For 



example, Zaheer and Zubairy have used path-integral methods |T9[, while several authors 
have used perturbative techniques The drawback of these latter methods is their 

restriction to small coupling. Finally, Swain has given formally exact result for the time 
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evolution of the Rabi Hamiltonian However, the solutions are expressed as continued 
fractions, which limits their usefulness in practice. 

In this paper we shall investigate the evolution of the Rabi system from an initial state 
composed of an empty field mode and an unexcited atom, which we denote as |0, |) and 
henceforth call the unexcited vacuum state. Within the RWA, the system would remain 
in this state indefinitely. However, inclusion of the counter-rotating terms means that we 
allow "energy non-conserving" processes to occur. It is these virtual processes that drive the 
evolution of the system. For simplicity, we shall only consider the resonant case, u = ujq. 

In order to investigate the evolution of the Rabi system from the state |0, |) we employ 
here one of the most versatile and most accurate semi-analytical formalisms of microscopic 
quantum many-body theory, namely the coupled cluster method (CCM) [23-30]. Coupled 
cluster techniques are widely regarded as being amongst the most powerful of all ab initio 
quantum many-body methods. As such they provide a number of distinct advantages over 
more traditional or more specialised methods which have hitherto been used in quantum 
optics and allied fields and, more specifically, to study Rabi systems. The CCM exists in 



two versions, the so-called normal (NCCM) and extended (ECCM) types |26|, and in the 
present paper we employ only the former version. 

A particular advantage of the CCM, a nonperturbative method originally developed in 
nuclear physics by Coester and Kiimmel [^, is that it has been extremely widely applied and 
tested on a huge variety of physical systems [^. These include areas as diverse as nuclear 
and subnuclear physics, quantum field theories (both in the spatial continuum and on the 
spatial lattice), condensed matter physics, quantum magnetism, and quantum chemistry. In 
almost all such cases CCM techniques now provide numerical results that are either the best 
or among the best available. 

A pertinent, but quite typical, example is the electron gas, one of the most intensely 
studied of all quantum many-body problems. Here, the CCM results p5|, |3l], |3^ for the 
correlation energy, for example, agree over the entire metallic density range to within less 
than one millihartree per electron (or < 1%) with the essentially exact Greens's function 
Monte Carlo results that are, very exceptionally, available for this fermionic system. The 
CCM results have never been bettered by any other technique. Elsewhere, in quantum 
chemistry, for example, after its early introduction by Cizek |]53[, the influence of the CCM 
has been profound, to the point where it it nowadays the method of choice for most highly 
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accurate chemical studies 11341 R5 



Since the Rabi Hamiltonian involves a two-level system which is modeled by Pauli pseu- 
dospin operators, we note finally that the CCM has also been very successfully applied 
to a large number of spin-lattice systems exhibiting anti-ferromagnetic and other forms of 
magnetic ordering. In instances such as unfrustrated models on bipartite lattices where a 
Marshall-Peierls sign rule [PS| exists (which provides a means to circumvent the infamous 
"minus-sign problem" inherent in simulating many-fermion systems by quantum Monte 
Carlo (QMC) techniques), the CCM provided results fully comparable in accuracy with 



those obtained by QMC means fS^]- The CCM is even able to predict with good accuracy 
the positions of the phase boundaries which mark the quantum phase transitions between 
states of different quantum order. Furthermore, the CCM provides equally accurate results 
for systems which are frustrated either geometrically (e.g., on a triangular lattice) or 
dynaimically (e.g., by competing interactions on different bonds on the lattice) where 
QMC methods are much more difficult to apply. In several such instances the CCM results 
are now the best available. 



Although the time-dependent formalism of the CCM |26| has existed for a considerable 
time, it has not yet found the same range of applications as its static counterpart. For 
example, Hoodbhoy and Negele have investigated the application of the technique in 
nuclear dynamics, using two interacting Lipkin Hamiltonian systems as a test model. It 
has also found more application in chemical physics [^, |T|. For example, Monkhorst 
has outlined the application of the time-dependent CCM to the treatment of molecular 
eigenstates with the aim to describe such phenomena as scattering, chemical dynamics, 
and laser chemistry. Sree Latha and Durga Prasad have investigated the application of 



the technique to multi-mode systems with vibronic coupling a mechanism describing, 
amongst other things, non-adiabatic phenomena in the electronic spectroscopy of polyatomic 
molecules. We note that all of these studies have tended to only use the lowest levels of 
CCM approximation, due to the complexity of the systems studied. In comparison, we shall 
use high levels of approximation in an endevour to describe the Rabi system as accurately 
as possible. 

For other applications of the CCM the interested reader is referred to the reviews con- 
tained in Refs. [^, An important feature of all the applications is that the results 
provided by the CCM are able to be systematically improved upon via well-defined hier- 
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achies of approximations. At each level of approximation the results are open to ready 
physical interpretation in terms of the correlated many-body clusters involved and their cor- 
responding multiconfigurational creation operators. To date the NCCM has already been 
used to investigate the ground-state of the Rabi Hamiltonian where it has been shown 



to give excellent results for values of the coupling parameter g < 0.665. In view of the 
demonstrable success in treating so many other many-body systems, we now wish to apply 
it to the dynamical evolution of Rabi systems, with a particular aim to shed light on the 
importance of the counter-rotating terms dropped in the commonly made RWA. 

II. FORMALISM 

We now briefly describe the application of the normal coupled cluster method to the Rabi 
Hamiltonian of interest to us here. 

Let 1^1' (t)) and (^E'(t)| be the exact ket and bra states at time t for an arbitrary many-body 
quantum mechanical system, chosen so that 

(^(t)|^(t)) = l, Vt. (6) 

The Hilbert space for our system may be described in terms of a model state or cyclic vector 
I'I'o) and a corresponding complete set of mutually commuting multiconfigurational creation 
operators Ic"]!. The set is defined with respect to the reference state, such that 

C/|$o) = = (^olC*], V / 7^ 0, in a notation in which Cq = 1, the identity operator. In 
general, J is a set index and the operators C\ contain products of single-particle operators. 
The set index {/} is complete in the sense that the set of states jcjl^o)! provides a complete 
basis for the Hilbert space. The reference state, |$o); must be chosen to be non-orthogonal 
to the actual wavefunction of the system, 

($o|vI/(t))^0, Vt. (7) 

Usually one chooses ($o|C/C'j|'^'o) = ^u^ where 5ij is some suitably defined Kronecker 
symbol. However, our later choice of l^*]! will leave the set |c]|<I>o)| orthogonal but not 
normalised; 

($o|C7C)|$o) = NiSij. (8) 
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The ket and bra states are formally parametrized independently in the normal coupled 
cluster method as 



where 



1^0 ly^O 



(9) 



(10) 



and k{t) is a c-number. Using this parametrization, the expectation value of an arbitrary 
operator, X, is given by 



{X) = {^it)\X\^{t)) = {%\Se-'Xe'\^o). 

In the evaluation of {X), we use the nested commutator relation, 

1 



(11) 



X = e-^Xe^ = X + 



x,s 



+ 



2! L 



x,s 



s 



+ ■ ■ 



(12) 



This equation describes a similarity transformation of the operator X. The similarity trans- 
formed Hamiltonian, H = e~^He^ , lies at the heart of the coupled cluster method. The 
ability of the NCCM to describe the Rabi stationary ground state has been investigated 
elsewhere |||, i5]. Here we shall only consider the time-dependent formalism. To this end, 
we introduce the action functional 



A 



' \m\(4.-H]\^{t))\. 



to 



The stationarity principle. 



dt 



6A ^ 5A 



(13) 



for all independent variations in the bra and ket states such that (5|\l/(tj)) = 
5(\l/(tj)| ; j = 0, 1, reproduces the time-dependent Schrodinger equations. 



(14) 



H\^) = ^-|-^\^), (vl/|i7=_,|(vl;|. 



(15) 



Equation (11^ gives Hamilton's equations of motion for the cluster coefficients. 



dt Ni dsi 
where Ni is the norm from Eq. (P). 



.dsj_ _ 1 d{H) _ .dsi _ 1 d{H) 



dt Ni dsi ' 



(16) 



Following Ref. |Q we make the following choice of reference state and operators, 

|$o) = |0,i), S = Si + S2, (17) 

with 

oo oo 



n=l n=l 



The corresponding expansion for S is 

§ = 1 + li + I2, (19) 

with 

00 00 
Si = Y.^rl\t)b\ t = Y.^:\tW-'a~. (20) 

n=l n=l 

We note that the states (fe^)" |0, |) and (t"!")^" cr^|0, i) are eigenstates of the excitation 
operator, with eigenvalues equal to n. They have corresponding even or odd parity (i.e. 
eigenvalues of 11 equal to +1 or —1, respectively) depending on whether n is even or odd. The 
reference state |$o) has positive parity. This is the same as the ground-state, as demanded 
by Eq.(l^). The overlaps, Ni, of the operators (|T^) are 

iVW = (0,i|6"(6t)"|0,i)=n! 
iVf = (0, 1 |a-6"-i(6l")«-V+|0, i) = 4(n - 1)! . (21) 

This choice of reference state and cluster operators has been found to predict a spurious 



phase transition aX g ^ ^ 0.665 |T0[. This in fact signals a breakdown of the calculation 
for g > gc- For g < g^ this scheme gives excellent results for the ground state, and so it is 
within this region that we shall work. 

III. TIME EVOLUTION 

We shall start the system in the unexcited vacuum, |0, |). This corresponds to the initial, 
t = to, conditions 

4'Hto) = s^\to) = = ^:\to) = ^\to). (22) 
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Other choices are certainly possible, but the above choice is the most obvious one that 
corresponds to a physical state and which satisfies Eq. (|^) . 

In the evaluation of the time evolution equations we truncate the sums in Eqs. ( ]T8| ) and 
(PPP at n = A^, giving the so-called SUB- approximation, in which all the coefficients s£' (t) 
and s^n\t), i = l,2JoTn>N are set to zero. The resulting set of coupled equations (|TB[) 
is then solved numerically by taking finite time-steps. We note that since the initial state 
|0, J.) has even parity 11, all the coefficients Sn\t) and s^\t), i = 1,2, vanish identically 
for odd values of n, at all times. The convergence of our method with decreasing step size 
is demonstrated in Table |. Naively, one would expect a similar type of convergence with 
increasing SUB-A^ level. Adding more cluster operators should allow the actual wavefunction 
to be described more accurately. However, we find that this is only the case for small 
coupling, as shown in Table ||. In general we find that for a given coupling parameter, 
g, there is a maximal value of the truncation index A^ for which convergence occurs. For 
higher values of A^ the method diverges, with the calculated cluster coefficients diverging 
to infinity. To understand this behaviour we look at the energy spectrum calculated by the 
NCCM under this scheme. This spectrum is determined by first calculating the ground-state 
coefficients, {sj, s}}, and using these to construct the so-called dynamic matrix, Hd of linear 
response theory . This matrix is then diagonalised to give the energy eigenvalues. A full 
discussion of this procedure would take us too far affeld and so we shall just consider the 
results and refer the reader to Ref. for further information. Fig. |T] shows the excitation 
energies as determined by a SUB-4 calculation. A SUB-A^ calculation yields AN eigenvalues, 
half of which are spurious negative energy solutions, generated as a result of the symmetry 

The important feature to note is that at certain values of g (e.g., at g ^ 0.35 for the 
SUB-4 case in Fig. |lD, which depend on the truncation index A^, two real energy levels come 
together, and as g increases, become a complex conjugate pair. Thus, for higher values of g, 
the NCCM predicts a pair of complex conjugate energy eigenvalues for the Hamiltonian. We 
can see what this means for the time evolution by considering the basic quantum mechanical 
expression for the evolution of an arbitrary wavefunction, 

1^(0) = 5^(w„|^(to))e-^^"(*-*«Vn), (23) 

n 

where \un) is the eigen-ket corresponding to energy eigenvalue En- We see that having 
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Step size 



0.0333 
0.0250 
0.0200 
0.0100 
0.0050 
0.0033 
0.0025 
0.0010 
0.0005 



-9.168663x10-3 
-9.167811x10-3 
-9.167578x10-3 
-9.167426x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 



-1.772499x10 
■1.772437x10 
■1.772420x10- 
■1.772409x10 
-1.772408x10 
-1.772408x10 
■1.772408x10 
■1.772408x10" 
■1.772408x10' 



(2) 



2.946230x10-^ 
2.365733x10-8 
2.212346x10-8 
2.114473x10-8 
2.108448x10-8 
2.108127x10-8 
2.108073x10-8 
2.108048x10-8 
2.108048x10-8 



^ < s 



,(2) 



-2.517358x10-8 
-2.493248x10-8 
-2.482074x10-8 
-2.473925x10-8 
-2.473400x10-8 
-2.473372x10-8 
-2.473367x10-8 
-2.473365x10-8 
-2.473365x10-8 



TABLE I: Convergence of the method with decreasing step-size. The values of the real and 
imaginary parts of two typical coefficients are shown at time gt = 1, for a SUB-12 calculation with 
g = 0.05. 



3? 



Ml 



-9.050635x10-3 
-9.167532x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 
-9.167416x10-3 



■1.783754x10" 
-1.772413x10 
-1.772408x10" 
■1.772408x10" 
■1.772408x10" 
■1.772408x10" 
-1.772408x10 
-1.772408x10" 
■1.772408x10" 
■1.772408x10" 
■1.772408x10" 



3? 



1.675505x10- 
2.126271x10" 
2.108048x10"" 
2.107987x10- 
2.107987x10" 
2.107987x10" 
2.107987x10" 
2.107987x10" 



4.392506x10-9 
-2.456155x10-8 
-2.473365x10-8 
-2.473357x10-8 
-2.473356x10-8 
-2.473356x10-8 
-2.473356x10-8 
-2.473356x10-8 



TABLE II: Convergence of the SUB- A'" method with increasing truncation index N. The values 
of the real and imaginary parts of two typical coefficients are shown at time gt = 1, for g = 0.05. 
The step size was 0.0005. This convergence is only valid for small couplings. 



10 



Figure 1 
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FIG. 1: The excitation-energy spectrum of the resonant Rabi Hamiltonian (cj = = 1) as a 
function of g, as determined by a SUB-4 NCCM calculation. Solid (dotted) lines show real (imag- 
inary) parts of the excitation energies. The negative-energy solutions and those with imaginary 
components are spurious. 

imaginary components to En of both signs will lead to exponentially growing terms in 
this sum, as opposed to bounded oscillations. Thus we expect our NCCM time evolution 
calculations to break down for a given coupling if the SUB-A^ level predicts complex energy 
eigenvalues at that coupling. Fig. ^ shows the values of g, for a given SUB-A^ truncation 
index A^, above which the spectrum contains imaginary components. These are the points 
where our time-evolution calculations break down. The SUB-2 time evolution breaks down 
at exactly the point where the ground-state calculation does. 

Sree Latha and Durga Prasad have encountered similar problems in their applications 
of the closely related multi-reference time-dependent coupled cluster method (MRTDCCM) 
P6| , which also posits an exponential Ansatz for the wavefunction, but uses a model space 
spanned by a number of states as opposed to the NCCMs single reference state. They trace 
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Figure 2 
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FIG. 2: The critical vahic gc = gdN) of the couphng parameter, above which the NCCM time- 
evolution calculations break down, is shown for various values of the SUB-iV truncation index A''. 
This is the coupling above which the SUB-TV excitation spectrum develops imaginary components. 
The Hamiltonian is resonant (a; = = 1), 

the origin of the complex eigenvalues back to the use of the similarity transformed Hamil- 
tonian H, which, when truncated, is not necessarily Hermitian. This lack of Hermiticity 
permits complex eigenvalues, whose corresponding eigenvectors are related to the the so- 
called "intruder states" . They conjecture that this situation may be expected to arise when 
the reference state or space interacts strongly with a state or set of states in the rest of the 
Hilbert space; which is just the situation we observe here. 

In understanding the following results, it will be useful to study a Fourier transforms of 
the time series. Terms similar to 

m = f2{k-i)\sf\t)^^\t) (24) 

k=l 

occur in all the quantities that we look at below. We define F{0,) as the Fourier transform 
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FIG. 3: Fourier transform of the quantity f(t) from Eq. ( ^41) for the resonant Hamiltonian 
{uj = loq = 1) with g = 0.05. This shows a discrete spectrum with three main frequency components. 



of f(t) and Fig. ^ shows a plot of this quantity for typical parameters. We see a discrete 
spectrum with three main peaks. These peaks correspond to the lowest three positive-parity 
energy levels in the spectrum, as one would expect. This means that the time evolution will 
be quasi-periodic. This quasi-periodic structure is reflected in the behaviour of the cluster 
coefficients, as can be clearly seen from the typical parametric plot shown in Fig. ^. 



IV. ATOMIC INVERSION 

The atomic inversion, (cr^), has been the primary atomic quantity of interest when study- 



ing the Rabi system, not least because it is experimentally determinable |47|. In our NCCM 
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Figure 4 



0.035 




-0.035 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

-0.06 -0.05 -0.04 -0.03 ,,.-0.02 -0.01 0.00 0.01 

(2) 

FIG. 4: A parametric "phase-space" plot of the coefficient (t), shown for a range of the time 
parameter t for which < gt < 120. {g = 0.05, SUB-30, uj = ujq = 1) 

scheme the atomic inversion is given by 

oo 

(0 = -l + 8 5^(n-l)!4^)5l^). (25) 

n=l 

Being an observable, {a^) should always be real but the cluster-coefficients Sn^ and 
i = 1,2, are, in general, complex. The truncation of the cluster operators leads to the exact 
hermiticity of the bra and ket states being broken. This in turn means that the atomic 
inversion calculated under the NCCM is not constrained to be real. Calculations show that 
S2^'s2^ is the dominant contribution to the sum in Eq. (|25| ) and reveal these two coefficients 
to be almost complex conjugate to one other. Subsequent terms in the summation conspire 
to reduce the size of the spurious imaginary part, this reduction becoming more perfect 
with increasing truncation. For small couplings we can almost completely eliminate this 
imaginary part, although the restriction on the maximum SUB-iV level permissible for higher 
couplings means that we cannot eliminate this component entirely as g becomes larger. The 
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FIG. 5: Time evolution of the atomic inversion, (cr^), for two different couplings: a) g = 0.05, 
SUB-30 and h) g = 0.2, SUB-14. The Hamiltonian is resonant, (oj = ujq = 1). 

complex part of the calculated (a^) is small and can be used as a rough measure of the error 
in the NCCM calculations. For the parameters used here, RMS (5^ {(cr^)}) ^ 2 x 10~^^, 
Max (5^ {(cT^)}) ^ 2 X IQ-" for g = 0.05 and, RMS (55 {{a')}) ^ 1 x 10-^ Max (5> {{a')}) ^ 
4 X 10~^ for g = 0.2. Figure ^ shows the evolution of the atomic inversion for two different 
couplings, g = 0.05 and g = 0.2. 
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V. FIELD OBSERVABLES 



A. Photon Number 

The most important operator associated with the field is b^b, the photon number operator. 
In terms of the NCCM coefficients, it has an expectation value n{t), 

CO OO 

n{t) = m = Y,nn\ s^^^^^ + 4 ^(n - l)(n - 1)! (26) 

n=l n=l 

Figure |^ shows the evolution of this quantity. The nature of this evolution is very similar to 
that of the atomic inversion, and we clearly observe the atom exchanging energy with the 
field. A time average of (6^6) provides an estimate of the mean number of virtual photons in 



the cavity at any given moment. For example, for g = 0.05, (feTfe) ~ 0.0063 and for g = 0.2, 



{Wb) ^ 0.12. 



B. Photon Anti-bunching 

The next field observable we shall study is the second-order correlation function, 

(2), , (&t(t)6t(t + r)b(t + r)b(t)) 

^ {bKmt))^ ■ ^ ^ 

This allows us to study whether the field exhibits anti-bunching. For this purpose, we only 
require g^'^\0), 

a^fn) - ^^'^'^^^ (28) 



and we define the convenient parameter ^ 

y = (b^b^bb) - (b^b)^. (29) 
A value of y < corresponds to g^^\0) < 1 and indicates an anti-bunched field. 
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FIG. 6: Time evolution of the field occupation number (n) = {b^b), using the same parameters as 
in Fig. 5. The Hamiltonian is resonant, {to = ujq = 1). 
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For the NCCM, (b^b^bb) is given by 

oo 

{b^b^bb) = s«5<J) 

n=l 
oo 

+ ^ r2m(n + m)! 4^)5(^)51^1^ 

n,m=l 
00 

+ 4^(n-l)(n-2)(n-l)!.(f)5<f) 

n=l 

00 

+ 8 J] n(m - l)(n + m - 1)! 4'^s(^)sS„. (30) 

n,m=l 

Figure |^ shows the quantity y for various couphngs. For small values of g we see that y never 
drops below zero and thus the field is never anti-bunched. For higher values of coupling, we 
see that y occasionally does fall very slightly below zero. This effect is extremely small and 
we believe that it is due to inaccuracies introduced by using small truncation levels. 



C. Squeezing 

By analogy to the position and momentum operators of the harmonic oscillator, we 
introduce the following quadrature operators for the electromagnetic field ||50| , 



Qi = - [6exp {iut) + 6^exp {—iut)'\ , 

Q2 = — \bexp (iut) — b^ exp {—iujt)~^ . (31) 

Their variances, given by 

(AQO' = I {{b^b) + 3? [{b') exp (2^u;t)] } - {3? [{b) exp + i 

{AQ^f = i {(6^6) - 3? [{b') exp {2iu;t)] } - [{b) exp {lut)]}' + ^, (32) 

satisfy the uncertainty relation, 

AQ1AQ2 > |. (33) 
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Figure 7 
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FIG. 7: Time evolution of the photon anti-bunching measure, y, using the same parameters as 
Fig. 5. Note that y < indicates anti-bunching. The Hamiltonian is resonant, {lo = loq = 1). 
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The field is said to be squeezed when either (AQi) or {AQ2) < j. The NCCM expressions 
for the remaining expectation values in Eq. (|3^) are easily calculated, and are given by 



n=3 

00 00 

+ ^"'^^ + m-2y. 4^)s«^|^_2 (1 - 6n,i5, 



m,l I 

n=l m=l 

00 

+ 4$:(n-i)!.(2)e^- 



-2 

n=3 

00 00 



+ 8 5^5^n(m-l)(n + m-3)! .^.^^^ig^^, (34) 

n=l m=2 

00 00 

(6) = .« + $:n! «l, + 4$:(n-l)! .i^)^!. (35) 

n=2 n=2 

Figure || shows the evolution of these variances from the vacuum. For g = 0.05, we see 
that the values of (AQj)^; i = 1,2 only just drop below j; (min(A(5j)^ = 0.2495 for 
g = 0.05). However, for greater couplings, we do see that squeezing is more significant; 
(min(AQ,)2 = 0.1351 for g = 0.2). 



VI. DISCUSSION 

We have demonstrated the ability of the NCCM to describe the time evolution of a 
simple but important quantum system. We have also outlined the limitations of the method. 
Despite which, we have been able to obtain a range of useful reults for the system. 

We have seen that the counter-rotating terms give rise to quite complex behaviour in the 
evolution from the unexcited vacuum. The field, although not anti-bunched, does exhibit 
squeezing, which becomes more pronounced with increased coupling. However, the small 
absolute magnitude of these effects and the limitations of the model, such as neglect of 
thermal photons, clouds the experimental significance of the results. 

In assessing the performance of the NCCM in describing this system, it is useful to 
compare results with those obtained by the configuration interaction (CI) method using the 
same basis. The CI method is equivalent to diagonalisation in a truncated set of basis states. 
If we use the same number of states in both CI and NCCM calculations, the two procedures 
are of approximately the same computational complexity. In Table |T| we compare atomic 
inversions of the Rabi system calculated by both methods for two couplings. We have chosen 
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FIG. 8: Time evolution of the squeezing variances. Squeezing occurs when (AQi)^ or (A(52)^ < \- 
a) g = 0.05, SUB-30, short time, h) g = 0.05, SUB-30, longer time, and c) g = 0.2, SUB-14. The 
Hamiltonian is resonant (uj = loq = 1). 

times for this comparison where the atomic inversion is near a maximum, and thus the state 
system at these times is as far from the reference state as possible. It should be noted that 
both sets of = 16 results from the CI method may be treated as exact to the precision of 



Table III, as they are converged with results for much larger A^. These results demonstrate 
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g = 0.05; 


gt = 1.25 


g = 0.20; 


gt = 1.26 


N 


CI 


NCCM 


CI 


NCCM 


2 


-0.999999 


-0.981944 


-0.936100 


-0.745812 


4 


-0.981757 


-0.981759 


-0.696200 


-0.696457 


6 


-0.981759 


-0.981759 


-0.692675 


-0.694392 


8 


-0.981759 


-0.981759 


-0.693094 


-0.694211 


10 


-0.981759 


-0.981759 


-0.693087 


-0.692888 


12 


-0.981759 


-0.981759 


-0.693087 


-0.692909 


14 


-0.981759 


-0.981759 


-0.693087 


-0.692926 


16 


-0.981759 


-0.981759 


-0.693087 


-0.693384 



TABLE III: Comparison of NCCM SUB-A^ results with results obtained from CI diagonalisation 
including the same number of basis states. The table lists the atomic inversion of the resonant 
Rabi Hamiltonian [uj = ujq = 1) for two different couplings. The times, given by gt, where the 
inversion was evaluated, were chosen so that the inversion was near a local maximum. 



several things about the performance of the NCCM in this system. In the region where 
the NCCM spectrum contains no complex energies, the NCCM describes the system better 
than the equivalent CI diagonalisation, especially for low (A^ = 2,4) truncation levels. This 
is due to the superior counting of independent excitations in the NCCM [Q. Conversely, 
the g = 0.2 results reflect the fact that the presence of the complex energies prevents the 
NCCM from converging, limiting the accuracy of the NCCM for higher couplings. 

The initial aim of applying the NCCM to the Rabi Hamiltonian was to produce an accu- 
rate microscopic description of the time evolution of the system across the whole coupling 
range. This has however not been completely realized due to the incursion of complex ener- 
gies in the NCCM spectrum. Similar problems occur in the ground-state description. The 
underlying reasons for this failure of the NCCM are not yet entirely clear to us, and will 
merit further study. However, it does seem likely that the existence of the Juddian points, at 
which level crossings occur, indicates the presence of a subtle symmetry whose preservation 
or breaking is not reflected in our simple choice of reference state. It should be noted that 



although van der Walt |45| has tried a number of different reference states and operator 
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selections, none of these has yet entirely solved this problem. 

Future work includes analysis of the Rabi Hamiltonian with the NCCM in the holomor- 
phic representation, and an extension of the above method to evolution from a coherent 
state. 
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